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- - - Abstract 

cni: 

I In this paper we consider the secure transmission in fast Rayleigh fading channels with full knowledge of 

' the main channel and only the statistics of the eavesdropper's channel state information at the transmitter. For 



(N 



C/3 



the multiple-input, single-output, single-antenna eavesdropper systems, we generalize Goel and Negi's celebrated 
artificial-noise (AN) assisted beamforming, which just selects the directions to transmit AN heuristically. Our scheme 
may inject AN to the direction of the message, which outperforms Goel and Negi's scheme where AN is only 
injected in the directions orthogonal to the main channel. The ergodic secrecy rate of the proposed AN scheme can 
be represented by a highly simplified power allocation problem. To attain it, we prove that the optimal transmission 

■ scheme for the message bearing signal is a beamformer, which is aligned to the direction of the legitimate channel. 

m" 

00 ■ After characterizing the optimal eigenvectors of the covariance matrices of signal and AN, we also provide the 

CSJ ■ necessary condition for transmitting AN in the main channel to be optimal. Since the resulting secrecy rate is 

o: 

■ a non-convex power allocation problem, we develop an algorithm to efficiently solve it. Simulation results show 
that our generalized AN scheme outperforms Goel and Negi's, especially when the quality of legitimate channel 

><■ 

\^ ■ is much worse than that of eavesdropper's. In particular, the regime with non-zero secrecy rate is enlarged, which 

^ [ 

can significantly improve the connectivity of the secure network when the proposed AN assisted beamforming is 
applied. 
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I. Introduction 

In a wiretap channel, a source node wishes to transmit confidential messages securely to a legitimate 
receiver and to keep the eavesdropper as ignorant of the message as possible. As a special case of the 
broadcast channels with confidential messages [[II, Wyner [[3 characterized the secrecy capacity of the 
discrete memoryless wiretap channel. The secrecy capacity is the largest rate communicated between the 
source and destination nodes with the eavesdropper knowing no information of the messages. Motivated 
by the demand of high data rate transmission and improving the connectivity of the network [[3l, the 
multiple antenna systems with security concern are considered by several authors. With full channel 
state information at the transmitter (CSIT), Shafiee and Ulukus [[U first proved the secrecy capacity of 
a Gaussian channel with two-input, two-output, single-antenna-eavesdropper. Then the authors of [[5l-|[7l 
extended the secrecy capacity to the Gaussian multiple-input multiple-output (MIMO), multiple-antenna- 
eavesdropper channel using different techniques. On the other hand, due to the characteristics of wireless 
channels, the impacts of fading channels on the secrecy transmission were considered in [[51, [[H with full 
CSIT. Considering practical issues such as the limited bandwidth of the feedback channels or the speed 
of the channel estimation at the receiver, the perfect CSIT may not be available. Therefore, several works 
considered the secrecy transmission with partial CSIT [[9l- [[T3l . In [[9l- [[Tn . the authors naively chose the 
directions of signal and AN without optimization and the resulting performance is suboptimal. In addition, 
they solved the power allocation via full search, which is inefficient. Furthermore, they did not prove the 
equality of the power constraint is hold (using all power is optimal). In |[T2l . a single antenna system 
is considered, thus the authors did not solve the beamformer and power allocation problems. Also, the 
authors did not prove the rate increases with increasing total power. In lfT3l . the authors did not consider 
the AN in the transmission, and thus their scheme is a special case of ours. Indeed, as shown in [[9l, [[111 . 
adding AN in transmission is crucial in increasing the secrecy rate in fading wiretap channels. Also under 
the case that the main channel is fully known at transmitter, the optimal direction for signals is not solved 
analytically in [[T3l . However, the secrecy capacities for channels with partial CSIT are known only for 
some limited cases, i.e., the transmitter has single antenna with block fading [[TOl and only the statistics 
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of both the main and eavesdropper's channels are known at the transmitter lfT4l . 

In this paper, we consider an important type of wiretap channels with partial CSIT, namely, the multiple- 
input single-output single-antenna-eavesdropper (MISOSE) fading wiretap channels. We assume that the 
main channel has a constant channel gain and the eavesdropper channel is fast faded, respectively. We 
also assume that the transmitter has perfect knowledge of the main channel and only the statistics of the 
eavesdropper channel. We adopt the artificial noise (AN) assisted secure beamforming as our transmission 
scheme, where the AN is used to disrupt the eavesdropper's reception [fTTI [|9l. Although the secrecy 
capacity of the considered channel is unknown, the performance of the AN-assisted beamforming has 
been shown to be capacity-achieving in the high signal to noise ratio (SNR) regime when the transmitter 
is equipped with a large number of antennas [[TTIl . However, in other operation regimes, the heuristically 
selected directions in [[TT]| BH to transmit AN may not be optimal, where the AN is restricted to be in 
the null space of the legitimate channel. This motivates our study on optimizing the AN assisted secure 
beamforming. Note that the assumption that the statistics of the eavesdropper's channel are known at 
transmitter was also used in ^ to design the power allocation between the signal and the AN (see flU 
(8)]). Thus our comparison to the method in [j9j in Section |V] is reasonable and fair. 

The main contribution of our paper is that we propose a general AN scheme, which outperforms BU. 
More specifically, the optimal AN may be full rank under some channel conditions rather than low rank, 
as restricted in |l9l. In addition, we provide a simplified power allocation problem to describe the ergodic 
secrecy rate, which highly reduces the complexity of solving the rate. To attain it, we characterize the 
optimal beamforming directions and the power allocation strategies for AN. We also provide the necessary 
condition for transmitting AN in the main channel to be optimal. After characterizing the eigenvectors 
of the covariance matrices of signal and AN, the resulting rate becomes a non-convex power allocation 
problem and we develop an algorithm to efficiently solve it. Simulation results confirm that the full-rank 
AN provides rate gains over |l9l, especially through the enlarged non-zero rate region. Note that the secure 
connectivity in a network is assured by the non-zero secrecy rate of the transmitter-receiver pairs |l3l. Thus 
our scheme is very useful for the large scale wireless network applications, which is an important type 
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of applications of the MISOSE wiretap channels ||3l . 

The rest of the paper is organized as follows. In Section HI] we introduce the considered system model. 
In Section Ull] an intuitive explanation of the rate gain from the proposed scheme is provided. We then 
develop our main result, i.e., the ergodic secrecy rate, via three steps. In this section we also provide the 
necessary condition to have a full rank optimal covariance matrix of AN. In Section |IVl we provide an 
iterative algorithm to solve the power allocation problem. In Section |V] we demonstrate the simulation 
results. Finally, Section |VI] concludes this paper. 

II. System model 

In this paper, lower and upper case bold alphabets denote vectors and matrices, respectively. The 
superscript (.)^ denotes the transpose complex conjugate. |A| and \a\ represent the determinant of the 
square matrix A and the absolute value of the scalar variable a, respectively. A diagonal matrix whose 
diagonal entries are a\...ak is denoted by diag{a\ ...Qk). The trace of A is denoted by tr(A). We define 
C{x) = log(l +x) and (x)^ = max{0, x}. A-*- is the null space of A. The mutual information between two 
random variables is denoted by /(;). I„ denotes the n by n identity matrix. A >- and A ^ denote that 
A is a positive definite and positive semi-definite matrix, respectively, a >- b denotes a majorizes b. 

We consider the MISOSE system as shown in Fig. [H where the transmitter (Alice) has nx antennas and 
the legitimate receiver (Bob) and the eavesdropper (Eve) each has single antenna. The received signals at 
Bob and Eve can be respectively represented as 

yk = h"xk + ni^k, (1) 
Zk = gk^k + n2,k, (2) 

where G C"^^^ is the transmit vector, k is the time index, h is the constant main channel vector, 
gk ~ CN{0,lnj) is the random eavesdropper's channel, and nj^ and n2^k circularly symmetric complex 
additive white Gaussian noises with variances one at Bob and Eve, respectively. In this system model, 
we assume that full CSI of the legitimate channel and only the statistics of Eve's channel are known at 
transmitter. Without loss of generality, in the following we omit the time index to simplify the notation. 
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The perfect secrecy and secrecy capacity are defined as follows. Consider a (2"^,n)-code with an 
encoder that maps the message w E W = {1,2, . . . ,2"^} into a length-n codeword, and a decoder at the 
legitimate receiver that maps the received sequence y" (the collections of y over code length n) from the 
MISOSE channels ([T]) to an estimated message w E iV. We then have the following definition of secrecy 
capacity. 

Definition 1 (Secrecy Capacity 4701/ ).- Perfect secrecy is achievable with rate R if, for any positive £ 
and z', there exists a sequence of {2"^,n)-codes and an integer no such that for any n> no 

/(w;z'\h",g")/n<£,and Pr(H' ^ w) < e', (3) 

where w is the secret message, z", h", and g" are the collections of z, h, and g over code length n, 
respectively. The secrecy capacity Cs is the supremum of all achievable secrecy rates. 
From Csiszar and Komer's argument [[T]|, we know that the general secrecy capacity can be represented 

by 

C= max I{u;y)-I{u;z\g). (4) 

p(x|u),p(u) 

However, for our considered CSIT setting, which is not full CSIT, the optimal p(x|u) and p{u) are still 
unknown. We propose to apply the linear channel prefixing and Gaussian signaling to f{x\u) as 

x = u + v, (5) 

where u ~ CN{0, Su) and v ~ CN{0, Sy) are independent vectors to convey the message and AN, respec- 
tively. In addition, the feasible channel input matrices of signal and AN belong to the set 

S={(Su,Sv) :tr(Su + Sv) <Pr,Su^O,Sv^O}. (6) 

Substituting ([T]), Q, and ([5]) into dH), we have the ergodic secrecy rate with generalized AN (GAN) as 

^""^sS^rH i+h«s.h j-'^hi i+8Hs.g )\) ■ ™ 

Note that we do not limit the covariance matrix Sy of the AN v to have any special structure besides 
the conventional one Thus our GAN scheme generalizes the AN in BH, which is only allowed to 



be transmitted in the null space of the main channel. On the contrary, our GAN can be transmitted in 
all possible directions. We then solve the ergodic secrecy rate optimization problem ^ for the proposed 
GAN beamforming (GAN-BF) scheme in the following sections. 

III. Optimization of the ergodic secrecy rate 

In this section, we identify the structure of the optimal solutions S* and S* for the GAN-BF optimization 
problem (|7]), where AN is not restricted in the null space of the main channel. By exploiting the optimal 
structure, we transform the complicated optimization problem over the covariance matrices ([7]) as a much 
simpler one in Theorem [T] In the following Theorem [H the optimized ergodic secrecy rate of the GAN-BF 
is merely characterized by the power allocations among the message bearing signal, AN in the direction 
of the main channel, and AN in the directions orthogonal to the main channel. 

Theorem 1: For the MISOSE fast fading wiretap channel with the perfect information of the legitimate 
channel h, and only the statistics of the eavesdropper's channel g ~ C//(0,I„j.) known at the transmitter, 
the optimization of the secrecy rate in (|7]) can be reduced to the following optimization problem 
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Rgan = max 

PU'PVi-PV2- 

Pu+Pvi+{'''t~1)Pv2=Pt 





|h| 




l + lll 


^W^Pv, 



log 1+ . luM.n -E 



v 



log 



1 + , 



+ 



(8) 

where Pu , Py^ , and, Pv2 are the powers of the signal, the AN in the main channel, and the AN in the null 
space of the main channel, respectively. G, = \gi\^ ~ EXP{1), which is the exponential distribution with 
mean equal to 1 , for i = 1 , 2, . . . , nj. 

Comparing (|7]) to ([8]) we can easily find that the optimization problem is vastly simplified from solving 
two matrices to three scalar variables. Note that we divide the proof of Theorem \T\ into three parts for 
the tractability and each part corresponds to Theorem 2, Lemma 3, and Lemma 4, respectively. Before 
proving ([8]), we introduce two important lemmas to proceed. 

Lemma 1: Given a diagonal matrix D = diag{d\, d2,---, dn) G C"^". Assume di > d2> ■ ■ ■ > d,i and U 
is unitary. Then U = [h/| |h| |, h^/| |h| |] and U = [h-^/l |h| |, h/| |h| |] maximizes and minimizes h^UDU^h, 



respectively. 

Proof: We can rewrite the maximization problem in the statement of the lemma as 

n n 

max ^ di\hi\^, s.t. ^ \hi\^ = ||h||^, (9) 

i=l i=\ 

where h = U^h, hi is the rth entry of h. Then it can be easily seen that \hi \ = ||h|| with |/z2| = |^3| = 
•■■ = \hn\ =0 can optimize dH). Therefore, it is clear that U = [h/| |h| |, h-*-/! |h| |]. The minimization part 
can be proved similarly. ■ 

Now, we identify the eigenvectors of the optimal S* and S* through the following lemma. 

Lemma 2: The optimal covariance matrices of the signal and AN S* and S* for ^ have the same 
eigenvectors as [h/| |h| |, h^/| |h| |]. 

Proof: Assume Su + Sv and Sy are eigen-decomposed as UDiU^ and VD2V^, respectively. First, we 
can reform (fT4l) as 



/ /l+h^^UDiU«h\ r /l+g«UDiU«g\ 

maxR = maxmaxi? = maxmax log - — . l,,.^ — E log l,^.^ 

Su,Sv Di,D2U,v D,,D2U,v V Vl+h^VDzV^hy [ ^Vl+g VDzV^gy 



s, 

Since g is isotropically distributed 



(10) 



E 



log 



1+g^UDiU^g 
l+g«VD2V^g 



E 



log 



l+g^D2g 



which is independent of U and V. Thus the inner optimization problem on the right hand side (RHS) of 
(fTOl) becomes 

(U , V ) = argmaxlog [j^^jj^^^ )■ dD 

Then from Lemma [Hwe know that U = nu[h/| |h| |, h-*-/] |h| |] and V = nv[h/| |h| |, h-^/l |h| |] can simulta- 
neously maximize and minimize the numerator and denominator, respectively, where Hu and Hy are the 
permutation matrices such that the eigenvector h/||h|| is in the direction of the maximum and minimum 
entries of Di and D2, respectively. Therefore, R is maximized. As a result, Su and Sy have the same 
eigenvectors. ■ 
We then introduce the interlacing theorem in Lemma [3] [fTSl p. 182] which will be used in proving 
beamforming is optimal (Theorem [2l). 
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Lemma 3 (Interlacing theorem): Let M G C"^" be a Hermitian matrix and let a G C" be a given vector. 
We then have 

(a) lkiM±asi")<lk+iiM)<lk+2iM±asi^), k=l,2,...,n-2, (12) 

(b) ?i^(M)<X;t+i(M±aa^)<?i^+2(M), /:=l,2,...,n-2, (13) 

where ?iyt(A) is the kth eigenvalue of A in ascending order. 
First, we identify the rank property of the optimal Su*. 

Theorem 2: For the MISOSE fast fading wiretap channel with the perfect information of the legitimate 
channel h, and only the statistics of the eavesdropper channel g ~ CA^(0,I„^) known at the transmitter, 
with the proposed GAN-BF, the optimal co variance matrix of signal for (|7]) is S* = yj^hh^. 

Proof: Since the secrecy rate optimization problem © is non-convex, we can use the Karush-Kuhn- 
Tucker (KKT) conditions to find the necessary conditions for the optimal solutions. We first transform 
(17]) into the following form to simplify the KKT conditions 

^^^^ = UT.s^°^ [ 1+h^s.h J - ^ 

Compared with (|7]), in (fT4l) . we place the maximum inside the operation (.)+. The equivalence of dV]) and 
(fT4l) comes from the fact that we can represent Rqan by range of the objective inside ()+ in (|7]) as the 
union of the sets of positive and negative rates and R , respectively, as Rgan = max(i?+lj^ )^ = 
max(i?+, which is max(i?+) when i?+ is a nonempty set and zero, otherwise. On the other hand, 

(max(7?+(Ji?^))+ is also max(7?+) when i?+ is a nonempty set and zero, otherwise. Thus we know ([7]) 
and (fT4l) are equivalent. Let ?i > 0, \j/u ^ 0? and \|/v ^ be the Lagrange multipliers of the three constraints 



log 



l+g^(Su + Sv)g 
1+g^Svg 



(14) 



in Q, respectively, the KKT conditions of (|7]) is 



01 =s; = A(s:,s;) 

02 =s; = A(s;,s;) 



hh^ 

+ E 



l+h«S*h 



gg 



H 



Vus; =si\\f„ = 0, 
\i/vS; =s;\|/v = 0, 
tr(s;+s;) <Pt, s;^o, s;^o, 



(15) 
(16) 
(17) 
(18) 
(19) 



where 



A(S;;,S;)^aa^ + M, 
hh^ 



M 



l+h«(S*+S*)h' 



E 



gg 



H 



(20) 
(21) 
(22) 



Li+g^(s;+s*)g 

and S* and S* are the optimal input covariance matrices of u and v, respectively. In the following we 
denote A(S*,Sv) by A* to simplify the notation. After left and right multiplying ([H]) by (S*)-^, with 
([n]), we have the relation A*(Su*)^ = (Su*)^A* = XiS^f, where I = ''[^s^ ' P • Then we can apply 
|fT3l Lemma 8] to ensure X > 0, if i? > 0. Since A* and (S*)^ commute, they have the same eigenvectors. 
Therefore, we have 

Aa*As* = As*Aa* = ?iAs*, (23) 

where A^* and Ag* are the eigenvalue matrices of A* and S*, respectively. Due to M in (|20l ) is a negative- 
definite matrix |fT3l Lemma4], from Lemma [3l we know that all eigenvalues of A* are smaller to zero 
except for the largest one. This can be explained as following. By using Lemma [3] and letting k = — 2 
in (fT3l) . we have ?i„j._i(A*) < ?i„j.(M). Note that M is a negative definite matrix, i.e., ?i„j.(M) < 0. So we 
have ?i/(A*) < for ? = 1,2, ... — 1. Since X, is positive, from (|23] ) we know that it must be the largest 
eigenvalue of A*, i.e. X = Xnj{A*). In order to make the equality Aa*As* = ?iAs* valid, the eigenvalues 
of S* corresponding to non-positive eigenvalues of A* must be all zeros. Therefore, we obtain that S* 
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has only one nonzero eigenvalue. So the covariance matrix of S* is rank one if R> 0. Then with Lemma 
|2l we conclude the proof. ■ 
In the following we prove an important property, that is, using all the power is optimal for the proposed 
AN scheme. 

Lemma 4: To maximize (|7]), the sum power constraint in Q is hold with equality. 

Proof: Similar to Theorem 2, the key observation here is that with the selection of eigenvectors of 
signal and AN in Lemma [21 the first term on the RHS of (flOl ) is independent of the power of AN in the 
null space of the legitimate channel. Thus to find Py. for i = 2,3, . . . .nr given Py and Py^, the objective 
function becomes 



min E 

Pv2---Pv„j- 



( 



log 



GiPu 



1 + — 

. l + Gi/'v,+ EGA, 

\ i=1 



(24) 



From (1241) it can be easily seen that given Pu and Py^ , the value of the objective function decreases with 
increasing Pj. Thus we may change the first inequality constraint in ^ as an equality one. ■ 

Based on Lemma [21 and [H we have the following property for AN. 

Lemma 5: For the optimization problem ([7]), the optimal covariance matrix of AN is 



riT — I \ ||h||2 
Proof: To proceed, we transform ([24|) as 



1 fnTPv^-Pr+Pu^y^H^^p^_p^_p^^^, 



max E 



nj 



log I 1 + GiiV: + £g,/V, - log 1 + Gi {Pu +Pvy) + £gA 

i=2 J V i=2 



max Ejj 



(25) 



where the equality comes from the conditional mean, f{x) = E [log (a + x)- log (b + x)] and we denote 
1+Gi/Vj, l+Gi{Pif+Py^), and i^GiPy. by a, b, and jc, respectively. If given Gi =^i, Vgi, the optimal 

1=2 

power allocation of f{x) is Py^=Py^ = --- = Py^^ , then for the problem on the left hand side (LHS) of (|25l) . 
this power allocation is also optimal. This is due to the fact that G,- is unknown at transmitter by whom can 
not be used to change the power allocation. Therefore, we want to prove that under £ Py. = Pj—Pu —Pvi 

i=2 



iij 



riT- 



(26) 
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Here we introduce some results from the stochastic ordering theory [fT6ll to prove the desired result. 

Definition 2: (W, p. 234] A function \|/ : [0,°°) — )■ R is completely monotone if for all x > and 
n = 0, 1,2, ■ ■ • , its derivative exists and (-l)"\[;W(;c) > 0. 

Definition 3: |T6i (5.A.1)] Let Bi and B2 be two nonnegative random variables such that E[e^^^i] > 
E[e^^^-], for all s > 0. Then Bi is said to be smaller than B2 in the Laplace transform order, denoted as 

Bl <LT B2- 

Lemma 6: |[T6l Th. 5.A.4] Let Bi and B2 be two nonnegative random variables. If Bi <lj B2 then 
E[/(5i)] < E[/(52)], where the first derivative \[/ of a differentiable function / on [0,0°) is completely 
monotone, provided that the expectations exist. 

To prove (l26l) . we let 5i = £ G,/V,, ^2 = E GiPy to invoke Lemma [6l where Py. denotes the optimal 

i=2 i=2 

value of Pvi- It can be easily verified that \(^(x), the first derivative of f{x), satisfies Definition |2l More 
specifically, the nth derivative of '\\f meets 

^^T^ > 0, if n is even. 



i J^^ + W^<^^ if n is odd, 
when x> 0, since by definition, b> a> when R> 0. Now from Lemma [6] and Definition [3l we know 
that to prove ^ is equivalent to proving E[e"''^i] > E[e^''^^] or log(E[e"^'^i]/E[e"'''^2]) > 0, V5 > 0. From 
[fTTl p.40], we know that 

iiT riT 

£iog(i+2P4.)-£ 

k=2 k=2 

To show the above is nonnegative, we resort to the majorization theory ifTSl . Note that £^L2log(l +2Pvt>^) 
is a Schur-concave function in {Pv2, ■ ■ ■ tPv„^), > 0, and by the definition of majorization 



= i log(l +2^4^) - i log(l +2iV,^). (28) 

V*^[^ jy A:=2 k=2 



in.. 



p. . Pt-Pu-Pv, Pt-Pu~Pv, Pt-Pu-Pv, ^ 



we know that the RHS of (l28l) is nonnegative, V* > 0. Then (|26l) is valid. From Lemma |2] and |5l we can 
conclude that 



h/||h||,h^/||h| 



diag 



Pt-Pu-Pv, Pt-Pu-Pv,\ 

PVi , z , ■ ■ ■ , z 

Y nj — I nj — I J 



h/||h|Lh^ 



(29) 
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Then with the expansion 




I 



we conclude the proof. 



After substituting the S* from Theorem |2] and S* from Lemma [5] into (|7]), we can get ([8]). Note that 
when the main channel is fast faded but perfectly known at transmitter, as [[T2|. the achievable secrecy 
rate for this setting can be easily obtained from results in Theorem \T\ 



Although we have simplified the optimization problem in ([7]) to ([8]), since ([8]) is a non-convex stochastic 
optimization problem, it is still difficult to analytically solve the optimal power allocation Py, Py^, and 
Pv2 in dl]). Thus in this section we propose an iterative power allocation algorithm summarized in Table 
I, which can find solutions almost the same as the brute-force search. However, the complexity of the 
proposed algorithm is much lower than the one based on brute-force search. More specifically, the brute 
force search requires searching on a plane for the three variables Pu, Py^ , and Py^, simultaneously. However, 
the proposed algorithm divide the search into two sub-problems which costs much less complexity. Before 
introducing the iterative algorithm, we first provide a necessary condition in Theorem ([3]) for the optimal 
covariance matrix S* of the GAN to be full rank. This condition will be useful to test the correctness of 
power allocation found in proposed algorithm. 



IV. The iterative algorithm for power allocations between signal and 



GENERALIZED ARTIFICIAL NOISE 



First define 




where Ek{x) is the En-function [fT9ll . 



Then we have the necessary condition in the following. 
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Theorem 3: The necessary condition for the power allocation {Py ,Py^,Py^) to be optimal for ([8]) is 

+ (^ + ^^AiFi{Pv^) +A2F2{Pv,) 





1 


1 + 


||h| 


\'Pu 


1 + 1 


M^Pv, 1 + 1 


|h|| 




'u + Pv,) 



P \ B P A 

+ (nT-l + ^] t-^PkiPv.J - ^BnrFnAPv,) - {Pv, + {riT - \)Pv2)-^^Fi{Pu + Py,] 
PViJ k^\^V2 Pv2 Pu+Pvi 



' a' f P \ P 

\i-UPv +i'vj - - 1 + -^^ t.^kUPv2) + -^B,^F,,{Py,) ^ 0, (30) 
+Pyi \ PviJ k^i Pvi 



_PvA 
Pu 



then 

/A, A„ \ "7—1 D a' a' 

(^FiiPv.) + ^F2{Pv,)]+L ^FkiPv2)-^^F,{Pu+Py,)-—^F2{Pu + Pv 



"r-1 c' 

- I ^FkiPy,) ^ , , V' (31) 





|h| 


\'Pu 


(1.1 


|h| 


\'Pv) 


(l + ||h||2(Pf;+/Vi)) 



where 



X-riT Py. 1 ^^'^-'^n-Tv 



^1 = T ir^iw-^^ ^2 = — — T, 



TlT —k) ^- 



i.^iY'^Pvi' (i_^y-'' (^i_^v'""'+' 



, _ X-UT Pu+Py, . _ 1 j-^j 

with the requirement Py^ > 0. 

Now we present the derivation for the proposed iterative algorithm. The key idea of the proposed 
algorithm is as following. To prevent the high complexity of simultaneously solving Py, Py^, and Py^, 
we try to divide the problem as smaller ones and we can simply use bisection method to solve them. 
More specifically, we start from the KKT conditions, by eliminating the Lagrange multipliers, we form 
two equations each has different variables to solve. Then iteratively solve these two equations, we can 
find the power allocation. With the Lagrange multipliers X > 0, > 0, ^ui > 0, and fi2 > 0, by the KKT 
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conditions of ([8]), we then have 
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Assume that Pj), P^^, and Py^ are all non-zeros. Combining (|33l) . (|34l) . (|36l) . and (|37l) we have 
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Similarly, combining (|33l) . (|351) . (I36l) . and (|38l) . and using the fact that 
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since the channel gain of each antenna is independent and identically distributed (i.i.d.), we have 
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Now for the zth iteration, with a given Py^, we can find new {Pu,Py^) such that f2{PuTPy^,Pv',) =0 
according to (l4T|) . We can set Pu = {Pj — Pv2 — Pvi)/i^T — 1) then f2{Pu,Pvi,Pv2) becomes a function 
with only one variable Pv2- We let the resulted Py^ as Py^^'^ ■ Then with a given Py^^\ we can numerically 
solve a new Py^ such that fi{Py^,Py^^'')=0 according to (|39| ). We let the resulted as Py^^^^ and the 
iterative algorithm follows. The bisection method can be used to perform the numerical search. 

Based on the concept described above, we explain each step in Table |T] in detail. First, numerically 
finding the tuple {Py^^Py^.Pu) which exactly meet the equality (|39l ) (or (|4TI) ) is very hard. Therefore we 
relax (|39|) and (|4TI) by inequalities 

\fi{Py,,Py,)\<t, and \f2{Pu,Py,,Py2)\<tu (42) 

respectively, where ei is a small constant. Once the values from the bisection search validate the above 
inequalities, they are treated as the solutions of these inequalities. Together with the iteration step described 
in the end of the previous paragraph, we obtain Step 2 and 3 in Table IB Second, relaxing equalities (|39l ) and 
(|4T)) to inequalities (|42)) make solutions obtained depend on £i and may not satisfy the KKT conditions. 
Also the expectations in functions f\ and /2 ( (|39| ) and (|4TI) ) are calculated numerically via generation 
of the channel realizations. Thus as in Step 4 of Table HI we use the analytical results in Theorem [3] to 
verify the correctness of the solutions. Finally, the initial values for the first iteration in Step 1 are as 
follows. Note that two initial values are needed for specifying the search region of the bisection method. 
For initializing Step 2, the two initial values for Pu are and Pj—Py^, such that the corresponding values 
of function will have opposite signs. And there exists at least one solution in the interval [0,Pt — Py^]- 
By the same reason, for initializing Step 3, the two initial values for Py^ are and Pj —Py^inT — 1). In 
the /th iteration, the search regions are [0,Pr - /'If] and %Pt - {nj - l)/'!^ ] for and /i, respectively. 

However, the bisection method may not always work for searching solutions for I/2I < £1 in Step 2 of 
Table E Note that for the initial value Pu = Pr-Py}, fiiPr - Py^Py^ .0) < given P^f . On the other 
hand, given Py^, there exist two cases for /2 at initial value Pu = 0: one is that f2{0,Py^ ,Py^) < as 
depicted in Figure |2] (a), and the other is f2{0,Py^ ,Py^) > as depicted in Figure |2](b). In the later case, 
the bisection method works. However, if the former case happens, the function values have the same sign. 
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and the bisection method does not work. To solve this problem, we can use the golden section method 
II20I . which is a technique for finding the maximum in the interval [0,Pt — Py^^], i.e., to numerically find 
Pu first such that given f2{Pu,Pv^ jPy^^ is positive. After that we can follow the step 2 in Table I 
to solve Pu in the interval [Pu,Pt-Pv^]. If the maximum of f2{Pu .Py] ^Py]) in the interval [Q^PT-Py^] 
is still negative, we know that there does not exist any Py in this interval such that fi^Pu .Pyl .Py}) = 
given Py \ In this case, we set Pu =0 3& the solution of f2{Pu .Py^^ -^Py^) = given Py \ From simulation 
results, according to the iterative algorithm in Table H the power Py\ Py^ , and Py^ will converge to the 
optimal solution P^, Py^, and Py^, respectively, which satisfy the KKT necessary conditions. 

Remark 1: Note that in Section |IV] we assume that Py, Py^, and Py^ are all non-zeros to eliminate the 
multipliers. For channel conditions under which low rank AN covariance matrix is optimal, the proposed 
algorithm may have Py^ converge to a value approximately zero. When this value is smaller than a 
predefined threshold £2, we claim that Pvi = is optimal. 

V. Simulation results 

In this section, we illustrate the performance gain of the proposed transmission scheme over Goel and 
Negi's scheme. We use a 2 by 1 by 1 channel as an example. Assume that the noise variances of Bob 
and Eve are normalized to 1. From ^ we know that the rate Rqan only depends on the norm of the 
main channel. Therefore, we use | |h| p = 0.05, 0. 1, and 0.2 to indicate different channel conditions in 
the simulation. For the statistics of the eavesdropper's channel, we set E[Gi] = E[G2] = 1. In Fig. [3l HI 
and [51 which correspond to ||h|p = 0.05, 0.1, and 0.2, respectively, we compare the rates of Goel and 
Negi's scheme to that of our proposed signaling with the generalized AN. The blue and black curves 
represent searching the optimal power allocations exhaustively and by the proposed iterative algorithm, 
respectively. In the iterative algorithm, we set the iteration number MAX IT as 20, MAXCheck as 5, and 
£1 =£2 = 10"^ From Fig.HIl and [51 we can easily see that the proposed generalized AN scheme indeed 
provides apparent rate gains over Goel and Negi's scheme in the moderate SNR regions. In addition, we 
can observe that the rate gains decrease with increasing ||h|p, which is consistent with the results in [[T2l|. 
We can also find that the value of Pt which provides the largest rate gain also decreases with increasing 



17 

||h|p. This is because AN in the signal direction provides much more rate gains when Bob's received 
SNR is relatively small compared to Eve's. Furthermore, the power allocations of the proposed iterative 
algorithm indeed converges to those by exhaustive search. In and Fig. |6]we show the convergence rate of 
the proposed algorithm under ||h|p = 0.1 with different Pj. It can be found that the proposed algorithm 
converges fast under different Pt, i.e., it costs at most 7 iterations to the final value, which verifies the 
complexity of solving the power allocation is much lower than the full search. 

As another example, we also illustrate the optimal power allocation among Pu, Py^, and Py^ under 
||h|p — 0.05 in Fig. Ul It can be easily seen that as the received SNR increases, the power allocated to 
Pvj decreases and the rate gain over Goel and Negi's scheme also decreases. 

VI. Conclusion 

In this paper we generalized Goel and Negi's artificial noise (AN) for fast fading secure transmission 
with full knowledge of the main channel and only the statistics of the eavesdropper's channel state 
information at the transmitter. Instead of transmitting AN in the null space of the legitimate channel, we 
considered injecting AN in all directions, including the direction for conveying the dedicated messages. 
Our main result provides a highly simplified power allocation problem to describe the ergodic secrecy rate. 
To attain it, we proved that for a multiple-input single-output single-antenna-eavesdropper system with 
the proposed AN injecting scheme, the optimal transmission scheme is a beamformer which is aligned to 
the direction of the legitimate channel. In addition, we provided the necessary condition for the optimal 
covariance matrix of AN to be full rank. After characterizing the optimal eigenvectors of the covariance 
matrices of signal and AN, we also developed an algorithm to efficiently solve the non-convex power 
allocation problem. Through simulations, we verified that the proposed scheme outperforms Goel and 
Negi's AN scheme under certain channel conditions, especially when the legitimate channel is poor. 

VII. Appendix 

Before proving Theorem [3l we first introduce the following lemma which will be used. 
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Lemma 7: Given Di ^ D2, 
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Proof: We first write the expectation in (|43] ) in the following integral, 
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(45) 



\+Pv,t \ + {Pu+Py^)tJ (l+Py^t)"^ 

for i = 2,3, ... ,nT, and from [13, Lemma 4], we know that the non-diagonal entries of both the first and 
second terms of Y in dV]) are zeros, then Yij = for i ^ j. Therefore, we know that Y is a diagonal 
matrix and each diagonal entry from (l44l) and (|45l) is larger than zero, which completes the proof. ■ 
We now provide the proof of Theorem [3] 
Proof: We first rearrange (fT6l ) as 

02-C-Xl„^ + \|/[ = O, 

where 



C=UYU^-cc^, 
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;i+h^S;^h)(l+h^(Su + S*)h) 
Similar to (l23l) . we have 

Ac As; = As; Ac = tr(CS;)As;. 



(46) 
(47) 

(48) 
(49) 



And we know that the necessary condition for the optimal AN to be full rank is that when tr(CS*) > 0, C 
does not have any negative eigenvalues; or, when tr(CS*) < 0, C does not have any positive eigenvalues. 
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To verify this property, we resort to the fact from |fT3l Lemma 5] that if all eigenvalues X of aa^ — A are 
negative, then /(O) > 0, where is defined as, 

/(X) = l-a^(A + ?iI„r)"^a, (50) 

and A >- 0. Note that l(k) is a strictly increasing function when X > 0. Note also that C in (|46l ) is negated 
of aa^ — A. Thus all eigenvalues of C are positive implies /(O) > 0. Thus by substituting c and UYU^ 
into a and A, respectively, we have 



By Lemma |7] we know (UYU^) ^ exists. Then we can expand /(O) > from (|5T| ) as 

1 >c^ (UYU^)"^c. 

Then after substituting c from (|48] ) to (|52l ), and using Theorem [2l and Lemma [2l we have 

(l + ||h|pPv,) (^l + \\h\\^{Pu+Pv,i 
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From |fT3l Lemma 4] we know that Y is diagonal. In addition, with Y is invertible from the proof of 
Lemma Ul we can further rearrange the above as 
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Then by the definition of Y in (|47]) , and the fractional expansion, we can further express the above as 
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(53) 



where Ai , A2, A j , Aj, 5yt, and B'j, for k = 1,2, . . . .nj — I are defined in the statement of the theorem. In 
addition, tr(CS*) > implies 
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After some arrangement, (|54l ) can be further represented by 

+ (l + ^) A,F,{Py,) +A2F2{Pv,) + (nr - 1 + ^) I^^^^^^) 

(p \ 'IT p 
nT-\ + ^] t^'kFkiPv,) + -^B,^^F„^{Py^) > 0. (55) 
PViJ k=\ 
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TABLE I 

The iterative algorithm for power allocation between signal and generalized an 
Step 1 Set i = 0, = 0, and initialize search region for the bisection method. 

Step 2 Given P^^' and the total power constraint ([6}, find /V, (and thus Py = {Pt - Py^- Pv} ) I {"t - 1)) 
such that |/2(/'[/,/'y|',/'v,)| < ei, where /2 is defined in ^A\\ . 
Setp(;+"=Pv, 

Step 3 Given -P^'^'' and the total power constraint (|6}, find /Vi 

such that |/i(/'V|,/'y'^'')l < ei, where /i is defined in l [39l > 

setp(;+"=/v,. 

Step 4 Let / = /+ 1 and repeat Step 2 to Step 3 until MAXIT. 

Step 5 Check the whether the final power allocations meet Theorem [3] 

If not, randomly re-initialize /y^' and run Step 1-4 until MAXCheck. 
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Fig. 1. System model. 
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3. Secrecy rate versus transmit power under ||h|p = 0.05. 





Fig. 5. Secrecy rate versus transmit power under HhH" = 0.2. 




Fig. 7. Power allocation among Pu, Pv,, and Py, under t|h||2 = 0.05. 



